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Abstract. A well-known result from Hardy and Ramanujan gives an asymp- 
totic expression for the number of possible ways to express an integer as the 
sum of smaller integers. In this vein, we consider the general partitioning 
problem of writing an integer n as a sum of summands from a given sequence 
A of non-decreasing integers. Under suitable assumptions on the sequence A, 
we obtain results using associated zeta-functions and saddle-point techniques. 
We also calculate higher moments of the sequence A as well as the expected 
number of summands. Applications are made to various sequences, including 
those of Barnes and Epstein types. These results are of potential interest in 
statistical mechanics in the context of Bose-Einstein condensation. 

1. Introduction 

The first significant ideas dealing with the theory of partitions can be attributed 
to Euler [16] and were published in 1748. The next important milestone of partition 
theory was laid in 1918, when Hardy and Ramanujan [17j . using quite involved com- 
binatorics, produced their celebrated theorem, which gives the asymptotic result 
for the number of ways to express an integer as the sum of lesser integers. 

In the mid-1970's, the area of asymptotic analysis blossomed with two of its pri- 
mary works, those of Dingle 6 and Olver [28 . Among the plethora of information 
in these two books are some nice results concerning the method of steepest descent, 
more commonly called the saddle-point method. Also at this time, what is now 
one of the primary texts of partition theory, Andrews', The Theory of Partitions 
[3], became available. In 1975-76, Richmond used ideas of generating functions [24] 
and saddle-point methods to describe the moments of certain types of partitions 
[29] [30]. In the first of this two part series, The Moments of Partitions 2, using 
asymptotic analysis instead of combinatorics, Richmond reproduced the results of 
Hardy and Ramanujan [17] as well as calculating higher moments and variance. In 
his second paper, The Moments of Partitions II, Richmond went on to give asymp- 
totic results for the moments of a general sequence whose associated zeta-function 
has only one singularity in the interval (0, 1]. 

Though Richmond's results are remarkable, he does not consider the situation 
where the sequence gives rise to a zeta-function with arbitrarily many singularities 
at arbitrary values. Building upon Richmond's results, and taking full advantage of 
the saddle-point method as well as asymptotic analysis, the present paper addresses 
that question. 

As literature on zeta-functions is readily available, we give only those definitions 
and propositions which are used. Background on the Ricmann zeta-function can be 
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found in [9] and [32] • Though the Barnes [U [5] and Epstein zeta- functions [HI [15] 
are less well-known, they are of paramount importance in mathematical physics; 
see, e.g., Actor Q][2], Dowker [TUB], Elizalde QUI DUES and Kirsten [H [23], for a 
full reckoning. Details of a general zeta- function are given, e.g., by Voros [33] . 

After presenting the background, we derive what we call the General Moment 
Theorems. The first of these theorems addresses the question: given a sequence A 
of non-decreasing integers, with limited assumptions, how many ways are there to 
write a large natural number n as a sum of members from the sequence A? The 
other General Moment Theorem addresses higher moments. We then apply the 
General Moment Theorems to find the expected number of summands. 

We make applications to a variety of sequences, starting with recreating the 
results of Richmond's first paper, among other things re-deriving the famous Hardy 
and Ramanujan theorem described above. The General Moment Theorems are then 
applied to a sequence whose associated zeta-function has only one singularity. In 
essence, this is what Richmond describes in his second paper, except we assume 
no restriction on the location of the singularity. We then reproduce the results of 
Nanda [25] for two and three-dimensional Barnes-type partitions, and proceed to 
consider two and three-dimensional Epstein-type partitions. 

The advent of mathematical physics, specifically the area of statistical mechanics 
in the context of Bose-Einstein condensation, has brought partition theory to the 
forefront of current research, see, e.g., [19l [20] [26] . It is because of this context 
that we require the sequence A to consist of non-decreasing integers instead of 
increasing integers. More specifically, we allow the sequence to contain the same 
integer a finite number of times. In the physical context mentioned the repetition 
of the same number represents different quantum mechanical states with the same 
energy, so-called degenerate states. Distributing a given amount of energy among 
the states corresponds exactly to the partitioning problem considered here, where 
the repeated numbers correspond to the different states with the same energy. We 
hope that the results of this article will be of interest to not only the physics 
community, but to number theorists as well. 



2.1. Basic zeta-functions. In this section, we give the definitions and some prop- 
erties of specific zeta-functions needed for use in later sections. 

Definition 1. Let s G C with Sfts > 1. The Riemann zeta-function is defined as 



It is well known that a meromorphic extension of Ciz( s ) to the whole complex 
plane can be constructed and that the only pole of (ii(s) is at s = 1 [S]. 

Definition 2. Let s G C with 5Rs > d, and c E R, r G R. d such that c + rhf > for 
all m G Nq- The Barnes zeta-function is defined as [H[5] 



2. Zeta-functions 




tiSN 




'0 



For c = it will be understood that the summation only runs over m G Nq 



d o~{0}. 
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For r = (1, 1, . . . , 1, 1) := 1 we will use the notation £g(s,c) := Cg(s, c|l), in 
which case we have the following expansion [H [5] : 

Proposition 1. Let r = l and c > 0, then 

oo 

(b(s,c) := Ce(s, C |l) = ^e[ d) (c + 

1=0 

where 

id) _ (l + d-l 
6 ' ~\ d-1 

In the case r = 1, the above result indicates that the Barnes zeta- function can 
be represented in terms of the so-called Hurwitz zeta-function. 

Definition 3. Let s £ C with 3?s > 1 and c e R with c > 0. The Hurwitz 
zeta-function is defined as 



(h(s,c) = V - — — - 
(n + c) s 



n=0 

The Hurwitz zeta-function can be meromorphically continued to the whole com- 
plex plane and its only pole is located at s = 1 . 

Barnes showed the following relation between the Hurwitz and Barnes zeta- 
function [H[S]: 

Corollary 1. The Barnes zeta-function Cs(s,c) can be meromorphically continued 
to the whole complex plane and for c € R with c > one has 

Cs(s,c) = f a^^ B^ k (c)C H (s + 1 - M- 



(fc - l)!(d- k)\ 

Here B\ a \c) := _B J - ct ' , (c|T') are £/ie generalized Bernoulli polynomials |27| defined by 

{-l) d ^ (-t) n - d 



Similar results can be derived for the case c = 0. For example in dimension d = 2 
one finds 

(i) Gj a (*,o) = Cw(«-i) + &(*), 

whereas in dimension d = 3 the answer reads 

Cb 3 M) = \ (( n (s - 2) + 3Ctj(s - 1) + 2Crc(s)) . 

Finally we will consider sequences related to sums of squares of integers; the fol- 
lowing zeta-function will be useful. 

Definition 4. Define Q(m,r) = r\m\ + r^m^ + ■ • • + r^rn^. Let s E C with 
s Rs> f , and c € R, r € R d such that c + Q(rh,r) > for all fh S Nq. The Epstein 
zeta-function [Ml [15] is defined as 

C£(s,c|f)= V , —Jr^ ^, s - 
Z~t. d (c + Q(m,r)) s 

In case c = 0, it is understood that the summation ranges over meNjj- {0} only. 
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For all cases meromorphic continuations of £s(s,c\r) can be constructed; see, 
e.g., [HI HH [15l [22] • Because of the complicated appearance, involving series over 
Bessel functions, we do not display them explicitly. 

2.2. A general zeta-function. In this section we first introduce the general type 
of sequences we want to use, giving them certain basic restrictions, then we examine 
residues and particular values of the zeta-function associated with the sequence. 

Throughout this paper we denote by A a sequence that is a nondecreasing se- 
quence of natural numbers such that the following hold. 

(i) 1 g A. 

(ii) The partition function, 



(2) 9(i) = 



AeA 

converges for t > 0. 
(iii) For t — ► + , Q(t) admits a full asymptotic expansion 

(3) e(t)~£^»*S 

neNo 

where i n G R with i n +i > i n , io < and where i n — > oo as n — > oo. Later 
on, we will occasionally also use the notation — i n = //„. 

The fact that we restrict A to be a sequence of natural numbers is due to the fact 
that we analyze partitioning problems of natural numbers. In writing down the 
restrictions (ii) and (iii) we follow [33], where it is shown that these requirements 
lead to the well-defined spectral functions considered in the following. In particular, 
with these restrictions on the sequence A of numbers, we are now in a position to 
define a general zeta-function (of A-type) as the following. 

Definition 5. Let A be a sequence as described above, and let s G C with 5fts > /zq. 
We define the general A-type zeta-function as 



ca(*) = Et 



Typically, zeta functions are built from eigenvalues of an elliptic (pseudo) differ- 
ential operator [TTJ [3T] , but it remains a perfectly viable spectral function in the 
given context [53] , 

For our later considerations, we will only need to know residues and particular 
values of (\. To find these, we will use the standard integral representation 

1 r°° 

Ca(s) = / t^e^dt. 



r(«) jo 

For the residues, and values of Ca( s ) at s = — n, only the small-i behavior of the 
integrand is relevant and we may focus our examination on the function 

r ( s ) Jo 

realizing that the residues, and values of Ca(s) at s = —n, are precisely those of 
L(s). 
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Formally substituting J3]) into the above equation and performing the integration 
yields the expression 

L(s) = TTT) f <S " £ ***** = TTT) £ A -^T 

This formal calculation can be made precise and so allows to conclude the following 
propositions [3T1 [33] : 



Proposition 2. For i n ^ k S No, the residues of Ca(s) occur at s = —i n , and 
furthermore, 

(4) Res {Ca(*)} = =^hr. 

»=-t„ l(-«nj 

Proposition 3. For n £ No, 

(A(-n) = (-l) n n\A n . 

3. General Moment Theorems 

3.1. General results for all k. We start this section by defining moments of 
partitions. 



Definition 6. Let pA{n,m) be the number of partitions of n into m summands 
where each summand is a member of A. The A;-th moment of p\{n, m) is denoted 
by t A (n), and is defined by 

where p\(n, m) = for m > n and p\{n, 0) = 0. 

For small n, the above definition is sufficient for finding values of t\(n) by hand 
or using a little computer program. For large n, the calculation is much more 
difficult; we must find a more reasonable way to compute t\(n). We resolve this 
problem by first constructing a generating function, then we continue by evaluating 
the coefficients of this generating function; see [2U HH] ■ 

Let us start by defining 

(5) G A (x, z) = J] (1 - z^)- 1 = £ PA(n, m)x n z m , 

AeA n6N mGN 

and 

Note that 

^ fe G A ( a ;,z)| z=1 :=^G A (x) = ]T ^ ^V(«,m)x" = ^ *i(n)x n , 

n£N meNo riGN 

so that we have constructed a generating function for t\{n). 

Since A contains only integers, only n € No occurs in the summation. We may 
therefore apply Cauchy's formula for Laurent series coefficients, so that for e > 
suitably chosen, we have 

*a(«) = ! # k G A (x)x-( n +Vdx, 

2m Jc{0,e) 
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where (7(0, e) is the circle of radius e about x = 0. With the substitution x = e a , 
this easily becomes 

(6) A{n) = ^-.[ e n ^ lo ^ kG ^ a »da, 

2m Js.p. 

where s.p. indicates a closed path that goes through the saddle-point a = a,}, of the 
integrand. 

The saddle-point a — is found as a solution to the equation 

(7) A( a+ I log ^ A(e - O )) =0 . 

As we will see in the following, the large-n expansion of the moments results from 
a small-|a| expansion of the saddle-point equation. 

In order to evaluate ((6]) as n — > oo, that is, as |et| — * 0, it will be necessary to find 
a more explicit form of the saddle-point equation. We first simplify d k G A {e~ a ), 
making two cases, namely k = 1 and k > 2. For k = 1, 



d 

dG A (x, z) = z — 
oz 



G A (x,z) — 



zx 

7 

AeA 

For convenience we denote 



Iaga J 

AeA / \AeA 

A 



(8) S A (x,z) = ]T and &Sk{x,z)\ s=1 :=&Sk(x). 

AeA 

From the above string of equalities, we have -dG A (e~ a ) = G A (e~ a )S A (e~ a ). 
Repeating the above process for k > 2 yields ^\ 

d k G A (e- a )^G A (e- a )S A k) (e- a ), 

with 



(9) S[ k \e- a )=j2 



6i!6 2 !---6 fc ! V 1! 



x • • ■ x 



2! J \ k\ 



where the summation is over all solutions b\, b%, . . . , € No of b\ + 2&2 + ■ • • + 
kbk = k (see [21] for more details on sums like ([9])). For readability, we denote 
the summation in © by ^(i?°S" A (e^ a ), t ? 1 S" A (e- Q ), . . . , ^ fc " 1 S , A (e~ a )). Using the 
above quantities, the saddle-point equation ([7|) reads 

f a+-logG A (e- a ) + -logS' A k) (e- a )) = 0, 



rfo 



n n 
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Figure 1. The contour from a — ioo to a + ioo in the complex plane. 



or, more explicitly, 
(10) 



e aX - 1 c( fe ) ( P -a\ 
AeA J A v e J 

We next show that for large n the solution to this saddle-point equation is unique. 

At this point, we begin to make systematic use of the following easily shown 
identity. 



z~ l T(t)dt, 



Proposition 4. Let a > 0, 6 > 0, and | argz| < \ — S. Then 

2m 

where the limits of integration define the contour shown in Figure QJ 
Applying Proposition 2] to the first term in (fTUj) . we have 



\ i pcr+ioo 

which gives 



AeA ' AgA ZgN 



1 



V— — = — ' 

^ e aA - 1 27ri 

AgA 



(T+ZOO 



a-\ A (t-lXn(t)r(t)dt, 



where the contour is similar to Figure [TJ but a > hq + 1 so that all of the residues of 
the integrand lie to the left of a. This range for er guarantees absolute convergence 
of the resulting series and allows for an interchange of summation and integration. 
In the right half-plane, the integrand has simple poles at t = /xq + 1, Mi + 1, ■ ■ • , 1, 0, 
and we find 

ai) y ^— = y CTC0 \t V ft + - + o {- 

V ' A^ e a\_l L^t a ^ + l ^ f« a \ a 

AgA i=0 v 

where the A^s are defined by ((3]). The term o(l/a) summarizes subleading contri- 
butions as a — > 0, the leading one of those behaving like l/a 1_e , e > 0, e depending 
on the location of the right most pole of Ca(s) on the negative real axis. 
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The analysis of the second term in (I10|l follows along the same lines. For k — 1, 
using again Proposition [4j we write 

(12) s A ( e - a ) = — / o-^AWCwWrWdt. 

For fc > 1 we need to focus on the terms "d* 3 " 1 S\(e~ a ) for s > 1. Note that ([5]) 
gives 

(is) r-^ A( ,) = £ E ^f'f = E4 S) E tt^w = E W*). 

AeA l± ,T ^ j=l ASA [1 X > j=l 

where S^ 3 (x) = J2\eA (i-x y ) s ' anc ^ tne c j S ' are defined3as in by = 0, — 
1, 4 2) = 1, 4 2) = 0, and for s > 2, 



( S +D = I jcf + (s-j + 1)4% l<j<s 
1 0, i = s + 1. 



Also, note the identity 

(14) £ c M = (s -l)!. 

Using as before the substitution x — e _a , we are interested in sums of the form 

AeA ^ ; 

A little arithmetic yields 

AeA v ' AeA /eN AeA/eN 

where e| s - ) = ( ' S ^ J . Applying Proposition [4] gives 



(T+ZOO 

a _t A _t (/+.7) _ *r(t)dt 



AeAieNo 

/"(T+ZOO 



\AeA / Ve 

Since here we have to deal with Barnes zeta-functions of different dimension, wc 
adopt the notation that (g is the Barnes zeta-function of dimension d. With this 
new notation we find, recalling Proposition [T] 

-| />(7-Hoc 

(15) S ^(e- a ) = — a-\ A (t)&(t,j)T(t)dt. 

^ 7r * Jo — ioo 

Although the small- |a| expansion of (fT2| and (fl5|) can be, and will be, obtained 
later on, at this stage let us content ourselves with the following observation. The 
leading a — > behavior of M (e~°), for Zfc > suitable, is seen to be of the form 



1 Correcting a typo in |29| 
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k , respectively a k log a, depending on the location of /io. We therefore will 



have 



(16) 



±S{ k) (e-) _ Cl (fc)+ C2 (fc)/logq 
(e- a ) a 



the numbers c\(k) and 02(h) depending on the location of fiQ. In all cases, it is 
seen that the fc-dependent correction to the saddle-point equation is of the order 
0(l/a). This will turn out to be of great importance for the analysis to follow. 

Using the small-|a| expansion displayed in (fTTj) and (|16[) . the saddle-point equa- 
tion (|10p allows us to uniquely determine a in terms of n, at least for large n. 

With the saddle-point known in terms of n, for large n, we determine asymp- 
totic answers for the moments by applying a theorem of Olver (Thm 7.1, p. 126 of 
[HI) to ©. Noting that 



jl(a+ilo g ^ A (e-°) 



1 dn 
n da 



this theorem takes the following form: 



Lemma 1. For ak the solution of (|10p . as n — > 00, we have 



(17) 



t\{n) 



e nak d k G A (e- ak ) 
2^ 




As indicated, has to be thought of as being replaced by its large-n asymptotic 
expansion so that (JTTJ) represents a large-n asymptotic expansion. This will be 
explicitly done once we start looking at specific sequences A in Section 4. 

The next important observation is that to leading order as n — > 00, the result 
for the moments is independent of the k used for the saddle-point ak- 

Proposition 5. For all k, as n — > oo, we have 

t k A (n)=tl(n)-S[ k) (e-"°)[l + (l)]. 

Proof. Applying equation ^7} amounts to evaluating e ™a+io g G A (e- a ) j (e~ a ), 
and dn/da at the saddle-points a = ak- We first note that the saddle-point equation 
(flT)]) . together with ((TTJ) and (fl6|). imply 



1 



ak \Cn(fJ-o + l)A_ w u 



where o(l) denotes terms that vanish as n — > 00. In particular, to leading order, 
S^\e~ ak ) and dn/da\ a=otk are independent of the saddle-point used. 

To show this independence for e nQ fe+ lo g G ( e k ) to the relevant order, consid- 
erably more work is necessary because of the exponential magnifying factor. In 
order to show the Proposition, we need to show that fiat + logGA(e~ Qfc ) = nan + 
logGA(e~"°) + o(l) such that the difference due to the saddle-point chosen only 
produces subleading order corrections. First, again from equations (fTU|) , (fTTj) and 
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If n \ ™ +1 



X 



(_ MO \ 
_l 1 ci(fc) + c 2 (fc)/ logn / n \ \ 

Mo + i Ctc(a*o + 1)^4-^0^0 vCtc(w> + i)^-mo^o/ "y 

where the leading fc-dependence of the saddle-point solution has been depicted ex- 
plicitly; ci(fc) is determined from Ci(fc) and c%(k). Using the product representation 
of G\(e~ a ) from we have 

— a\l i ro+ico 

(18) logG A (e- Q ) = = a -*a(*)fR(< + l)r(i)<ft- 

Depicting only terms relevant for the leading fc-dcpcndcnce of the contributions 
considered, we have 

IogG A ( e -°*) = a^Cnivo + l)A- m + o(a^°). 
Therefore, from equations (TO)) and (jTTJ) , we conclude 

na k + logG A (e- ak ) = ci(fe) + ^^c 2 (fc) + a^°(/x + 1)Ck(mo + l)A-^ +... 
= + ..., 

that is, there is no fc-dependence of this expression up to the order o(l), which 
shows the assertion. □ 

So in the following, saddle-point will always refer to a := a>o, which is a solution 

of 

AeA 

To exploit equations (JTTJ) and ([5]) we first need a more complete expansion of 
equation (fT8|) . The relevant integrand has simple poles at t = fa for i = 0, 1, . . . , d, 
and a double pole at t = 0. Calculation of the integral gives, as a — > 0, 

log G A (e~ a ) = £ ^L±^ A- M< - Ac log a + &(0) + o(l). 

i=0 a 

Since e" Q G A (e- Q ) = e »«+i°gG A (e- a ) ; we have the following lemma. 
Lemma 2. For a the solution of (Tl9|) . 



e na G A (e- Q ) = a- Ao exp 



,i=0 



[1+0(1)]. 



Changing focus to 4~ | =a, we note that this quantity is not exponentiated in the 
solution of t\{n). Thus, as mentioned, to obtain asymptotic results we need only 
determine the leading order. With this in mind, from equation (|11|). as a — ► 0, 

dn 

da 



MMO + 1 Wo + l)^ n [1 + 0(1)], 
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SO that 



' -2tt 

dn I 

da \ a=a 



= (27r M o(Mo + 1)Ck(mo + 1)A-^)' 



x a^ +1 [l + o(l)]. 



With the main calculations behind us, we now begin to consider specific k values, 
making two cases; k = and k > 1. The proceeding calculation brings us to our 
first general moment theorem. 

3.2. General A-type: k — 0. To evaluate the 0-th moment of the general A- 
type, t\(n), we apply (fTT)) . with k — 0, along with Lemma[5]and (|2"D|) , to give the 
following theorem. 

Theorem 1. For A-type partitions, 

t° A (n) = (27r/i (/xo + 1)Ctj(mo + l)^-^) - ' a ^ +1 -+> 

Ms: + 1 



x exp 



C7i(/ii + l)A_ M< +4o + CA(0) 



[1+0(1)], 



where a is the solution of and i/ie 's are defined by ([3]) . 



Once a is replaced by its large-n asymptotic expansion, this determines the 
large- n asymptotic expansion of i A (n). 

Let us recall that i A (n) gives the asymptotic result for the number of partitions 
of an integer n over the sequence A. 



3.3. General A-type: Case k > 1. To evaluate the k-ih moment of the general 
A-type for k > 1, t\{n), note that we need only evaluate S^\e~ a ) and then apply 
Proposition [3 

It is evident from (J9j) that within the calculation for general k > 1, we will need 
the specific calculation for ■d°Sty(e~ a ) = S^{e^ a ). 
Using Proposition |U 



a-*CA(t)Cw(t)r(t)dt. 



To evaluate 5A(e _a ) we must distinguish between a few cases. The above integral 
has different values depending on whether /io < 1, (J-o — 1, or jM) > 1. We treat 
these cases independently. 

For fj, < 1, the leading pole is the simple pole of (n{t) at t = 1, and so 



(21) 



.s\(, ") = ^ [i + o^ 1 -^; 



o 



For /io = 1, the integrand has a double pole at t = 1, and 

(22) 5 A (e-«) = i(FP{C A (t)}-A_ 1 loga) [l + O^ 1 ^)], 

a V*=i / 

where FP{£a(£)} denotes the finite part (or constant term) of the expansion of 

£a(£) around t = 1. 
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For no > 1, the integrand has a simple pole at t = /io, and 
(23) S A (e- a ) = fofcoM-w. \i + o( q^w-i^o-wA 

Having completed the evaluation of S\(e~ a ), we now focus on the terms i? s_1 S , A(e _a ) 
for s > 1 . As shown previously, see (fl~3f and (IT5|) , the relevant quantity to consider 
is 

S S A j (e- a ) = 7T- / a- t a(*)Ci, a) (*,j)r(*)dt. 

Note that we need k values for s, namely s = 1, 2, . . . , k. Previously in this section 
we considered s = 1, now we turn our attention to the following three cases for 
s > 2: s < fio, s = /io, and s > /io. 

Let us remark, that although for each given case better error terms could be 
given, in the generality considered there are many cases necessary. Therefore, in 
order to make everything more readable, we refrain from doing so. 

For s < /io, the leading term comes from the simple pole of Ca(^) at t = /io, so 
that 

s s^ (e - a)= $H»o,3)A-» [1 + o(1)h 
Now we have, using the above and (fT^j) . 

s 

^- l S A (e- a ) = a-^'A^Ys^&i^j) [I + o(l)] . 

i=i 

For s = no, the leading term comes from the double pole of (\(t)(g\t, j) at 
t = /io- From (fTJ) it is immediate that 



E S {^(*,i)} = ^ 



and so 



S A^ a ) = ^ (fP{CaW} + A_ s Fp{c^(^)} 



+^y(lK»)-loga)J [l + o(l)]. 

Thus, with ((Til), 

r- 1 5 A (e- Q ) = ^ ((* - 1)! • FP {&(*)} + A_ s - log a) + 

+E c ? ) F = p{4 s) (*.i)M 

Finally, for s > /io, the leading term comes from the simple pole of Q; S (t,j) at 
i = s, so that 

^'( e ^)=a- s CA( S ) [1 + 0(1)], 

and so 

(24) r- 1 5 A (e-«) = (s - 1)! • a-*C A (s) [1 + o(l)] , 

where again (QJJ has been used. 
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We are now in a position to evaluate (e~ a ) as denned in ([9]); for the meaning 
of the notation S followed by a /c-component expression consult the paragraph below 
@. Again we must consider three cases: Mo < 1, Mo = 1, and Mo > 1. 

If Mo < 1, (EH) and (HU) give 

Si k \e- a ) = £ (Ca(1), Ca(2), 2Ca(3), . . . , (fc - 1)!Ca(*0) [1 + o(l)]. 
If Mo = 1, we have directly from |22|) and (12411 that 

S[ h \e- a ) = a- k (gP {Ca(*)} - A_! log a, Ca(2), 2Ca(3), ...,(*- 1)!Ca(*0) 

Finally, if Mo > 1, then the leading term of iS A (e _a ) comes from b\ = k. This 
implies 

S^(e- a ) = (S A (e- a )) k [l + o(l)], 

so that |23|) gives 

(e~ Q ) = a-^ k A k _, [CnM] k [1 + o(l)}. 
We may now apply Proposition [5] to give the following theorem. 

Theorem 2. for k > 1, and a the solution of (lll[) . i/ie following results hold. 

(i) //mo > 1) then 

t h A (n) = t° A (n) ■ a-»° k A% [{n(»o)] k [1 + o(l)}. 

(ii) If Ho = 1; fien 

t\{n) = i°(n) ■ oT fc £ (FP {Ca(«)} - A_! log a, 

CA(2),2C A (3),...,(fc~l)!C A (fc))[l + (l)]. 

(hi) // Mo < I? then 

t\{n) = f°(n) • a- fc £ (Ca(1), Ca(2), 2£ A (3), ...,(*- 1)!Ca(*)) [1 + o(l)]. 

We call Theorems [T] and [2l the General Moment Theorems. 
For k = 1, the above theorem yields the following corollary. 

Corollary 2. For A-type partitions, with a the solution of (JTTJ) and t\(n) as given 
in Theorem^ we have: 

(i) For mo < 1, 

t\(n)=t° A (n) -^-[1 + 0(1)}. 
a 

(ii) For Mo = 1, 

t\(n) = t° (n) • - (FP {Ca(«)} - A_x log a) [1 + o(l)]. 

a \t=i / 

(iii) For mo > 1, 

*aW = *aW qA10 [1 + o(l)j. 
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3.4. Expected number of summands. In this section we give expressions for 
the expected number of summands. 

Definition 7. The expected number of summands of a A-type partition of an 
integer n, denoted by m A {n), is defined as 

WA(n) = ^R- 

Lemma 3. The expected number of summands of a A-type partition of an integer 
n is 

m A (n) = S A (e- a )[l + o(l)}, 
where a is the solution of pip . 

Proof. This result is an immediate consequence of Proposition [5] □ 

In light of Lemma [3l a direct corollary of equations (|2T|) . (|22|) . and ([23]) . is the 
following result. 

Theorem 3. For a the solution of the following assertions hold. 

(i) If fio < 1, then 

m A (n) = ^ [l + o(l)]. 
a 

(ii) If Ho = 1; then 

m A (n) = - (fP{Ca(*)} - A-i log a) [1 + o(l)] . 
a \t=i I 

(iii) If Ho > 1, then 

^A(n) = aMo [1 + o(l)] . 



4. Applications 

In this Section we apply the General Moment Theorems of the previous section to 
a variety of special cases. We first reproduce the results of Hardy and Ramanujan. 
Then we will examine the case of one singularity at [i > 0. Furthermore, we examine 
multidimensional applications in Barnes and Epstein type sequences, where in each 
case we calculate higher moments and the expected number of summands. The 
results frequently use the notation E followed by a /c-component expression. We 
remind the reader that this notation is explained below (|9j). 

4.1. Hardy and Ramanujan. Let A = N such that Ca( s ) — Cn( s )- This is the 
case with one singularity at [i = 1 and we find the very familiar results of Hardy 
and Ramanujan [T7] for the asymptotic number of ways to express an integer n 
as the sum of lesser integers. The following theorems follow immediately from the 
General Moment Theorems, see also |29j . 

Lemma 4. For sufficiently large n, 

7T 1 3 

a = - — + 0(n"3). 
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Theorem 4 (Hardy and Ramanujan [17j). The asymptotic number of ways to 
partition n over N is 

t (n) = — ^e'V 1 ? [l + C(n"5) 



4V3r, 



Corollary 3. Let k = 1, then 
tl („) = -^e- 
Corollary 4. For k > 2, 



3 7 + lo S 



6n 



t fc (n) = 



1 



'6n 



^ 7 + log 



7T 

6« 



7T 



C^(2), 2C^(3), . . . , (fc - 1)!Ctc(A:) l + 0(n"i) 



Corollary 5. TTie expected number of summands of a Riemann type partition of 
an integer n is 



m(n) 



(in 



7 + log 



6n 



1 + in-? 



4.2. Zeta-functions with one singularity. Let us suppose that we have a non- 
decreasing sequence of natural numbers A, whose corresponding partition function 
Q(t) admits the full asymptotic expansion 



i=0 

in which ko < 0, and ki > for all i > 0. Define fi = —kg. Then the zeta- function 
associated with the sequence A has only one singularity on the positive real axis, 
namely at t = [i. Having met the assumptions of the General Moment Theorems, 
we proceed to apply them to the above sequence A. 

Note that within the previous section we could not easily solve for n in terms of 
the saddle-point a. Since we now have only one singularity, the leading orders can 
be determined easily and (fTTp gives 

i 

" -+o(l/n) 



C^(m + i)M- 



//. + 1 n 



The General Moment Theorems now yield the following corollaries. 
Corollary 6. For sequences A as described, we have 

_1 l-2/tp 2Ap-2-n 

t° A (n) = (2tt(^ + 1)) 2 (Ktj(m + 2( " +1) n x 



x exp 



n \ M+1 



(» + l)((n(n + l)A-»)—i +Ca(0) 



For the higher moments, according to Theorem [21 we need to distinguish three 
cases. 



Corollary 7. For sequences A as described, we have for k > 1 £/ie following. 
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(i) For fi > 1, 

, / 71 \ ^ t + 1 k i 

(ii) For n=l, 

-k 

x E (U - 4r log > Ca(2), 2Ca(3), • . . , (A - l)!CA(fc)) 

x [l + o(l)] 

(iii) For fi < 1, 

*a(«) = *a(«) ' n** 1 + 

x ]T (Ca(1), Ca(2), 2Ca(3), . . . , (k - l)!0v(fc)) [1 + o(l)]. 

Applying Theorem [3] gives the following corollary. 

Corollary 8. for sequences A as described, the following hold. 

(i) If fj, < 1, then 

i 

mA(n)=fa(i) Ub.+iM-J 11 + 0(1)1 ■ 

(ii) If \x = 1, i/ien 

(iii) If fj, > 1, then 

u 

4.3. Barnes type moments. We will now consider two special cases of Barnes 
type moments, the two-dimensional and three-dimensional cases, with r = l and 
c = 0. These results correspond to two and three-dimensional oscillator assemblies 
considered by Nanda pZ5] , 

4.3.1. Two-dimensional. For the two-dimensional case, (1111) gives 

V n / V3(2C TC (3))*y n§ 36n V 7 
We use (P) to evaluate relevant residues and 

& a (o,o) = -ilog27r + &(-i). 

We may now apply the General Moment Theorems to produce the following 
corollaries. 
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Corollary 9. For the two-dimensional Barnes type sequence with c — and f=l, 



4» = 



2tt 



x exp 



o At J 

2§ 



Cw(2) 



2*(Cw(3))* 



r" 3 



(Ctz(2)) 2 
12Cw(3) 



+ &(-!) 



l + O n 



Making the substitution n' = n/ (2^(3)), for ease of comparison, the above 
corollary corresponds to that of Nanda [35] Eq. (34)]. 

Corollary 10. For the two-dimensional Barnes type sequence with c = 0, f — 1, 
and k > 1, 

, 31-24 fc 



2tt 



x exp 



3(C^(3))3 
2f 



n 3 



V n 

Cr(2) 



2i(fo(3))* 12Ck(3) 



1 + O (n~3 



Note that the case fc = of Corollary |TD] is precisely tg (n) of Corollary In 
general this is true when the largest singularity is greater than 1; that is, in this 
paper, fiQ > 1. Thus when considering the three-dimensional Barnes case, we will 
give only the general result for k € No- But first we give the expected number of 
summands as a consequence of Theorem [31 

Corollary 11. For the two-dimensional Barnes type partitions with c — and 
r = l, 



mB 2 {n) = Ctc(2) 



1 + O (rC* 



4.3.2. Three-dimensional. For the three-dimensional case, we evaluate to give 



a = (3C*(4))*n _ * + 



n 2 



3*fo(3) 
,4(Ck(4))^ 

'8C tc (2)C K (4)-3(C k (3)) 2 



3432(C K (4))f 
Using (|2.1[) we evaluate relevant residues and 



32n 







Ilog27r+|c^(-l) + ^CK(-2)- 



We apply the General Moment Theorems to derive the following corollaries. 
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Corollary 12. For the three-dimensional Barnes type sequence with c = 0, f = 1, 
and k € No, 



*B S (") = 



(3fr(4))-* 

4<7T 

x exp 



(C*(3)) 
4Ck(4) 

Ck(2)- 



3C^(4) 

3(C^(3)) 2 
8C^(4) 



3Ck(3) 



3Ck(4) 



3C^(4) 

4 



1 + [n~* 



where O - S ( U (4)y> h 2 ( >7jl- i J + 2^n\ l )- 



4Ctc(4) 



Setting fc = 0, and making the substitution n" = nj (3^(4)), the above corollary 
corresponds to that of Nanda [25] Eq. (51)]. 
Now, Theorem [3] yields the following corollary. 

Corollary 13. For the three-dimensional Barnes type partitions with c = and 
r = l, 



™>B 3 {n) = Ck(3) 



3Ck(4) 



4.4. Epstein type moments. We now consider the sequence A = {n^ + n|- 
n 2 d \n 6 Nq — {0}}. Using the resummation formula [TS] 



E 

l= — oc 



1 CXJ 



;=-c 



it is easily seen that the associated partition function has the small-i asymptotic 
expansion 

d 



where the j4_« are given by 



1 



n = 1, 



With the help of ^ this determines the relevant residues. Additional quantities 
needed in the corollaries stated below can be extracted from [T2J, E2 O [TH [22] . 

We now approach two and three-dimensional Epstein type moments in the same 
way that we considered the Barnes type. 

4.4.1. Two-dimensional. Note that for the two-dimensional case, (jTTJ) gives 



( **C*(j) 
U>/2(Cw(2))* 



n * \-0 (n * ") 

8n V / 



Using this information we obtain the following corollaries to the General Moment 
Theorems. 



GENERAL MOMENT THEOREMS FOR NON-DISTINCT UNRESTRICTED PARTITIONS 19 



Corollary 14. For the two-dimensional Epstein type sequence with c = and 

f=T, 



4W = K( R (2))' 



An 



x exp 



^ + %^l- + 4(o) 



32Ck(2) 



1 + O n~ 



Corollary 15. For i/ie two-dimensional Epstein type sequence with c = 0, r = 1, 
and A; > 1, 



4(n) = (A B (2))" 



An 



x exp 



CL(o) 



E fP{& 2 W}-glog 



4?? 



, Cfi 9 (2), 2C& (3),..., (*-!)!&,(*;) 



x 1 + (n 3 



Theorem |3] gives the expected number of summands. 



Corollary 16. For the two-dimensional Epstein type sequence with c = and 
r = l, 



me 2 {n) 



An 



FP {fc(t)} _I log (E|ffl 



1 + n"* 



.<R(2) 

4.4.2. Three-dimensional. For the three-dimensional Epstein case, Ijlip gives 



'37riCw^ 




3*Ck(2) 




2hoo(c^(|)) ; 



5C.(f)c(|)-2(C(2))^-|-^0(,-f) 



Here we may note that /io = § > 1, so that we have the following corollary to 
the General Moment Theorems. 
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Corollary 17. For the three-dimensional Epstein type sequence with c = 0, r = 1, 
and k £ No, 




x exp 



O- - 

3 5 7T- r 




2*20 (C^(|))^ 

(C^(2)) 3 3C k (|)Ctj(2) 



n : > 



4(o) 



1 + (rT^ 



io(C«(f)) 2 °m- 

Theorem [3] yields the following corollary. 

Corollary 18. For the three-dimensional Epstein type sequence with c — and 
r = l, 

3 

tt'Ctc (I) / 16n 



mg 3 (n) 



5. Conclusion 



In this article we have developed a systematic approach to obtain, for large n, the 
moments of the number of partitions of n into m summands from a given sequence 
A of natural numbers. Under suitable assumptions on A, see ([31), the moments are 
expressed in terms of information obtained from the small-t asymptotic expansion 
of the associated partition function. The main technical tool is the saddle-point 
method in combination with the application of ([¥]). 

Whenever the small-t asymptotic expansion of the partition function can be 
easily determined, the large-rt expansions of the moments can be easily found. This 
has been demonstrated in Section 2] with sequences given by linear and quadratic 
sums of integers. 
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